Polaronic Quasiparticles in a Strongly Correlated Electron Band 
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We show that a strongly renormalized band of polaronic quasiparticle excitations is induced at the 
Fermi level of an interacting many-electron system on increasing the coupling of the electrons to local 
phonons. We give results for the local density of states at zero temperature both for the electrons and 
phonons. The polaronic quasiparticles satisfy Luttinger's theorem for all regimes considered, and 
their dispersion shows a kink similar to that observed experimentally in copper oxides. We calculate 
the quasiparticle weight factor z and deduce the local effective inter-quasiparticle interaction U. 
Our calculations are based on the dynamical mean field theory and the numerical renormalization 
group for the hole-doped Holstein-Hubbard model and large on-site repulsion. 
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I. INTRODUCTION 

Polarons as quasiparticles, in which the electronic ex- 
citations are strongly coupled to the lattice modes, have 
been studied theoretically for well over 50 years. Nearly 
all these studies have been confined to models with only 
one or two electrons, and do not include the electron spin. 
Such models have limited applicability and there is a need 
for studies of more realistic many-electron situations 1 . 
There are technical problems, however, in handling such 
models as non-perturbative methods are required to de- 
scribe the effects of strong electron-lattice interactions 
and few such methods are applicable to the many- 
electron case. Progress has been made in recent years 
in dealing with models with strong local interactions via 
the dynamical mean field theory (DMFT) 2 ^. Calcula- 
tions based on the Holstein and Holstein-Hubbard mod- 
els, where the electrons are coupled to local Einstein lat- 
tice modes, have been made using the DMFT to examine 
the effect of strong electron-phonon interactions on the 
Mott transition 4 ' 5 , for example. The Holstein models, 
however, was originally introduced to describe small po- 
laron effects. Within DMFT, only the single electron case 
can be solved exactly^. There have been few studies at 
finite electron densities and including the electron spin 9 , 
except in one dimension 10,11 . The problem is that if spin 
is included in the Holstein model, there is an effective 
attractive local interaction between the electrons, which 
leads to local bipolaron formation and, away from half 
filling, superconductivity. 

Once the electron spin is included the more realistic 
model is the Holstein-Hubbard model which includes an 
local electron repulsion U. In most situations the repul- 
sion U dominates and undermines any tendency to the 
formation of local bipolarons. In the Holstein-Hubbard 
model at half filling with a large U, the electron-phonon 
interaction plays no important role as charge fluctuations 
are suppressed; its main effect is to renormalise and re- 
duce the effective value of At half filling and small 
U, where the electron-phonon coupling dominates, again 
there is little evidence 4 of purely polaronic effects, be- 
cause of the formation of bipolarons 4,12 . The question 



arises, therefore, whether there is any parameter regime 
where purely polaronic effects can be observed. If there 
is, then does the behaviour correspond to the predictions 
based on the spinless models with only one or two elec- 
trons, or do many-electron correlations play a significant 
role? 

In this letter we present calculations to show clearly 
that small strongly renormalized polaronic excitations do 
exist in the Holstein-Hubbard model and form a distinct 
narrow band at the Fermi level. The Hamiltonian for this 
model is 
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where e(k) describes the dispersion of the band electrons 
and U is the on-site inter-electron repulsion as in the 
Hubbard model 13 . The electron density rii = n,f at 
site i couples linearly to the local displacement operator 
Xi = (b i + b\) / \J2mujQ with an electron-phonon coupling 
g. The phonons are assumed to be dispersionless (local 
Einstein phonons) with energy and m is the mass of 
the vibrating ions. 

Our calculations are based on DMFT using a semi- 
elliptical density of states po( w ) of the non- interacting 
system given by 
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(2) 



for \lu + po| < D and zero elsewhere, corresponding 
to a Bethe lattice, where 2D is the bandwidth and po 
is the chemical potential. The effective impurity prob- 
lem is solved using the numerical renormalization group 
(NRG) 1 ^ 5 ", as generalized for the calculation of dynami- 
cal response functionsiSiiLiS, and incorporating the den- 
sity matrix generalizationi 9 - due to Hofstetter. 

To suppress any tendency for bipolaron formation we 
take a large value of the on-site repulsion, U = 6 in units 
where D — 2, which for half filling corresponds to a Mott 
insulator^S. We choose a value of /i corresponding to 
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FIG. 1: Local density of states for the quarter-filled Hubbard 
model (U = 6.0) 

quarter filling (rij) = 0.5. The local density of states 
for g = is shown in Fig. ^ with the Fermi level such 
that the lower Hubbard band is 1/3 filled. In this regime 
the strong correlation effects are largely suppressed due 
to the charge fluctuations; there is no sharp quasiparticlc 
peak at the Fermi level, which occurs for the small doping 
regime near half filling, and there is a moderately reduced 
quasiparticle renormalization factor z rs 0.65. The small 
spectral weight between the upper and lower bands is due 
to the tails in the broadening in the higher energy peaks 
of the discrete NRG results. 

To study the impact of electrons coupling to phonons 
we concentrate on the lower Hubbard band. In Fig.[2]wc 
show the effect of increasing the electron-phonon inter- 
action g, taking luq = 0.2 in all cases. As g increases to 
the value g = 0.65 a sharp feature develops at the Fermi 
level which we attribute to a polaronic renormalization. 
Due to the large U, there is negligible double occupancy 
(riiinii) < 0.01 for this range of g indicating no local 
bipolaron formation. 

To look at the band aspects of these excitations, we 
introduce the k-resolved one-electron Green's function 
Gk,cr(w) with a local self-energy Yl a (u>), 



Figure shows a plot of the spectral density of the 
Green's function Gk(u>) for —2 < e(k) < 2 and g = 0.60. 
When integrated over k, it gives the local density of states 
for the corresponding g shown in Fig. [3 We plot the posi- 
tions of the maxima in the quasiparticle peaks as a func- 
tion of e(k) in Fig. Q] The flat dispersion in the region 
of the Fermi level is clear evidence of a narrow polaronic 
band (full line). We see a kink at u> ~ —0.05 which we 
discuss later in relation with the phonon spectra. None 
of this is seen in the pure Hubbard case, shown as the 
dashed line in Fig. 

The locus of points kp on the Fermi surface of the 
interacting system is given by 

e(k P ) = ii - 53(0) . (4) 
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FIG. 2: The development of the polaronic quasiparticle band 
in the lower Hubbard band of the quarter-filled Holstein- 
Hubbard model (U = 6.0) for increasing electron-phonon cou- 
pling g. 




-3-2-101234 



to 

FIG. 3: Spectral density pk(u>) of the electron Green's func- 
tion for the quarter-filled Holstein-Hubbard model (U = 
6.0,s = 0.60). 



According to Luttinger's theorem^! the volume of the 
Fermi surface should be the same as that of the non- 
interacting system. For a local self energy 53(w), this 
implies that the Fermi surface given by Eq. should be 
the same as the non-interacting one, e(kp) = fio- If these 
are to coincide, then /iq = fi — 13(0). 
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FIG. 4: Positions of the maxima of the spectral function pk{oj) 
of the electron Green's function for the quarter-filled Holstein- 
Hubbard model (U = 6.0, g = 0.60, full line) and the quarter- 
filled Hubbard model {U = 6.0,9 = 0.00, dashed line). The 
arrow indicates e(kp). 



To check the theorem we have calculated the density 
(rii) directly from the NRG expectation value in the 
ground state, and then from the non-interacting Fermi 
surface using the value of /io = (J. — S(0) from equation 
Q. The two sets of results are compared for a range of 
hole dopings in Fig. (left) for the pure Hubbard model, 
and in Fig. [S] (right) for the Holstein-Hubbard model with 
g = 0.5. The agreement between the two sets of results 
provide a convincing evidence of the applicability of Lut- 
tinger's theorem in strongly correlated systems both with 
and without the electron-phonon interaction. 

If we assume Fermi liquid theorySi and expand the 
self-energy in Eq. J2J about the Fermi level w = 0, and 
retain the first two terms, we can define a quasiparticle 
Green's function 



Gk,a(w) = 



1 



A-?(k) 



(•5) 



where e(k) = ze(k), and z is the usual wavefunction 
renormalization factor z = (1 — S'(0)) _1 , and jl = 
S(0)) is a renormalized chemical potential. The cor- 
responding density of states po{oS) for the non- interacting 
quasiparticles is given by 



irD 2 



D 2 -(lu + ,i) 2 



(6) 



where D — zD plays the role of a renormalized band 
width. We interpret this as a renormalized band of free 
polaronic quasiparticles. Integrating this free quasiparti- 
cle band up to the Fermi level gives the quasiparticle oc- 
cupation (hi). The resulting (hi) — (rii), demonstrating 
the one-to-one correspondence of the electron and quasi- 
particle states of Fermi liquid theor} ^ 22 ! 23 (see Fig- EJ). 

In Fig. we give the results for the change in the quasi- 
particle renormalization factor z with increase of g for 
quarter filling and U = 6. The values of z we calcu- 
lated by differentiating the self-energy, and also from an 
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FIG. 5: Density as a function of the chemical potential fi 
calculated directly (full lines) and from the volume of the 
Fermi surface (circles) for the Hubbard model (left-hand plot, 
U = 6) and the Holstein-Hubbard model (right-hand plot, 
C/ = 6, 3 = 0.5). 
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FIG. 6: Quasiparticle weight as a function of g for the quarter- 
filled Holstein-Hubbard model (U = 6.0) deduced from the 
derivative of the self energy (o) and from an analysis of the 
renormalized parameters (x). 



analysis of the renormalized parameters 2 ^^ on the ap- 
proach to the low energy fixed point. The two types of 
estimates are in reasonable agreement; the deviation for 
smaller values of g is largely due to fluctuations in the 
numerical differentiation of in this range. The be- 
haviour of z with increase of g is quite different from 
that near half filling where z increases with increase of 
g. Near half filling the strong correlation effects stem 
from the enhancement of spin fluctuations and the sup- 
pression of charge fluctuations due to the large value of 
U . In this regime the increase of z with g can be under- 
stood as due to the reduction in the effective value of U, 
UcS — U — 2g 2 hence an increase in g reduces U e g 
and hence reduces the renormalization^. The results 
for z at quarter filling shown in Fig.[|Jl by contrast, show 
a large decrease of z with increase of g. This increased 
renormalization is a pure polaronic effect; the coupling 
to the lattice causing an increase in the effective mass of 
the conduction electrons and a reduced mobility. 

The renormalized parameter approach^ 2 ^ also allows 
us to calculate the local interaction between the quasi- 
particles U. The results for U over this range are shown 
in Fig. There is a progressive reduction with in- 
crease of g as might be expected for the induced attac- 
tive interaction due to phonon exchange. For the largest 
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FIG. 7: Effective local interaction U between quasiparticles as 
a function of g for the quarter-filled Holstein-Hubbard model 
(17 = 6.0). 
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FIG. 8: Spectral density of the phonon propagator for the 
quarter-filled Holstein-Hubbard model (U = 6.0) for various 
values of the electron-phonon coupling. 



cupancy < 0.01) and their effective attraction is due to 
retardation effects. 

For values of g > 0.65 the dynamical mean field theory 
no longer converges satisfactorily, the iterations fluctuate 
between occupied and unoccupied states. This indicates 
the on-set of a broken symmetry state, probably some 
form of charge ordering or phase separation. This in- 
stability is also revealed in the one-phonon propagator 
d(u>) — ((^b i ;&!)) whose spectral density is shown in 
Fig. |21 With increasing g there is a progressive soften- 
ing of the phonon mode from its bare value luq = 0.2. 
For g w 0.65 considerable spectral weight develops for 
u> < 0. This weight relates to the expected number of 
excited phonons in the ground state. It clearly indicates 
an incipient lattice instability for g > 0.65. 

The kink to be seen at w = —0.05 in the quasiparticle 
dispersion in Fig. 0] is clearly related to the frequency of 
the renormalized phonon excitation at the correspond- 
ing g — 0.6 (Fig. ISJ . An analysis for different values 
of g shows a clear correlation of the kink energy to the 
renormalized phonon frequency. This is similar to the 
experimentally observed kink in the hole-doped copper 
oxidesSi. 

To summarize: We have shown that an accurate treat- 
ment of the Holstein-Hubbard model at finite electron 
density and T = reveals a clear polaronic quasiparticle 
band induced at the Fermi level. Luttinger's theorem is 
shown to be satisfied. A kink found in the quasiparticle 
dispersion relation is related to the renormalized phonon 
frequency. 



value of g the value of U becomes negative, even though 
U — 2g 2 /u>o > 0. This can be interpreted as due to the 
fact the very low energy quasiparticles hardly ever oc- 
cupy the same site due to correlation effects (double oc- 
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